Vortex nanoliquid in high-temperature superconductors 
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Using a differential magneto-optical technique to visualize flow of transport currents, we reveal 
a new derealization line within the reversible vortex liquid region in the presence of a low density 
of columnar defects. This line separates a homogeneous vortex liquid, in which all the vortices 
are delocalized, from a heterogeneous "nanoliquid" phase, in which interconnected nanodroplets of 
vortex liquid are caged in the pores of a solid skeleton formed by vortices pinned on columnar defects. 
The nanoliquid phase displays high correlation along the columnar defects but no transverse critical 
current. 

PACS numbers: 74.25.Qt, 74.25.Op, 74.25.Sv, 74.72.Hs 



Columnar defects (CDs) created in high-temperature 
superconductors by high-energy heavy ion irradiation 
act as very efficient pinning centers for vortices aligned 
along the CDs. The effects of this correlated disorder 
on the vortex properties have been investigated exten- 
sively mainly in the limit of relatively high density of CDs 
where the enhancement in the critical current for practi- 
cal applications is most significant and the vortex matter 
is known to form the Bose glass phase [1-4]. In con- 
trast, our understanding of the opposite limit, in which 
sparse CDs perturb a higher concentration of vortices, is 
still incomplete. Recent experimental studies [5,6] have 
shown that in this case the vortex matter can no longer 
be regarded as a homogeneously pinned medium. In- 
stead, an intrinsically heterogeneous structure with two 
distinct sub-systems of vortices with very different char- 
acteristic energies is formed: vortices residing on CDs 
are strongly pinned and create a rigid skeleton, while the 
interstitial vortices form relatively ordered nanocrystals 
intercalated within the voids of this porous skeleton. This 
experimental finding of a "porous" vortex solid has con- 
tributed to a series of recent theoretical studies and nu- 
merical simulations [7-12]. Moreover, it was found that 
the melting line B m (T) of the "porous" vortex solid has 
an unconventional behavior showing a pronounced kink 
[5,6]. This has led to the suggestion of the possible ex- 
istence of two different liquid phases: at fields below the 
kink a usual homogeneous vortex liquid is formed upon 
melting, whereas at fields above the kink the "porous" 
solid melts into a heterogeneous state in which liquid 
vortices coexist with a skeleton of pinned vortices. It 
was proposed that in such a case a new transition line 
separating these two liquid phases should be expected 
[5,3,10]. Using a new imaging technique we present here 
the first experimental evidence for such a transition line 



within the reversible vortex liquid phase that separates a 
heterogeneous "nanoliquid" from a homogeneous liquid. 
One of the interesting properties of this novel nanoliq- 
uid phase is the absence of in-plane dissipationless su- 
perconductivity, while the coherence along the c-axis is 
apparently still preserved. 
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FIG. 1. (color online), (a) DMO image of melting in a 
BSCCO crystal at B = 45 G, T = 79.75 K, and T modu- 
lation of 0.3 K. All the pristine regions of the sample are in 
the liquid phase (blue) while the circular regions irradiated 
with B<t> = 60 G are still solid (brown), (b) AFM 20 x 20 
/im 2 image of an etched mica irradiated along with BSCCO 
crystals (B<j> — 20 G). The dark spots reveal the locations of 
CDs while the bright yellow regions are examples of two voids 
in the distribution of columns. 

The CDs were created along the c-axis in 
Bi 2 Sr 2 CaCu208 (BSCCO) single crystals (T c ~ 90 K) 
using low irradiation doses of 1 GeV Pb ions at GANIL, 
France. The samples were covered with 50 /am thick 
stainless steel mask with a triangular array of 90 /am 
apertures. The mask blocks the ions and as a result CDs 
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are formed in the BSCCO crystals only within the 90 
/im apertures. Figure la shows the melting process in 
one of the BSCCO crystals of about 1.9 x 0.8 x 0.025 
mm 3 visualized by differential magneto-optical imaging 
(DMO) for magnetic fields H || c-axis. The first-order 
melting transition (FOT) is detected by observing a step 
in the equilibrium local magnetization [13-15], which 
appears as bright paramagnetic regions in the DMO im- 
ages [5,16]. For clarity, the liquid and solid regions are 
colored in blue and brown respectively. In Fig. la all 
the pristine regions of the sample have melted already 
(blue) while the 90 jam irradiated regions are still solid 
(brown). Figure 2 shows the B — T phase diagram with 
the corresponding melting lines in the pristine, £?^(T), 
and in the irradiated, £?^ D (T), regions with = 60 G 
(B^ = n0o 5 n is the areal density of CDs, and 0o is the 
flux quantum). The use of the irradiation mask allows 
accurate determination of the difference in the melting 
temperature in adjacent irradiated and pristine regions 
in the same sample. Figure 2 shows that the irradiated 
regions melt at substantially higher temperatures while 
strikingly preserving the first-order nature of the transi- 
tion [5] . With decreasing temperature along the B^ D (T) 
line the strength of the FOT in the irradiated regions 
diminishes, becoming weakly first-order and eventually 
turning into a second order transition at a temperature 
that is irradiation dose and sample dependent [5,17]. 
The open circles in Fig. 2 were obtained using field mod- 
ulation DMO while the solid circles using temperature 
modulation. 
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FIG. 2. (color online). B — T phase diagram in BSCCO: 
B^{T) - pristine melting, B^ D (T) - melting line with CDs 
with B^ — 60 G, and Bdi(T) — derealization line. The black 
lines are guides to the eye and the red Bdi(T) line is a fit 
to exp(— T/To). The arrow indicates the location of the kink 
in the B^ D (T) line. Inset: phase diagram of a B^ = 20 G 
irradiated sample. 

We now focus our attention only on the irradiated re- 



gions which are macroscopic on the scale of the inter- 
vortex distance. Figure lb shows an AFM image of an 
etched mica irradiated along with the crystals, which 
reveals the typical random distribution of CDs follow- 
ing Poisson statistics. The CDs form a dense network 
or matrix that contains numerous voids that are free of 
columns. Two such voids, or pores, are outlined schemat- 
ically by the light color in Fig. lb. The CDs act as strong 
pinning centers for the vortices and therefore at low B all 
the vortices reside on CDs. When the field approaches 
the matching field B^ most of the CDs become occupied 
by vortices. This network of pinned vortices forms a rigid 
skeleton, while upon further increase of £>, the extra vor- 
tices create ordered nanocrystals within the pores of the 
skeleton, as revealed by magnetic decoration studies [6]. 




FIG. 3. (color online), (a) DMO image of magnetic field 
self-induced by a transport current of 50 mA in the central 
part of the sample at B = 145 G and T — 74.2 K. The 
B^ = 60 G irradiated regions are in the nanoliquid phase. 
The bright and dark shades correspond to positive and nega- 
tive values of Btrans- (b) Profiles of Btrans along the yellow 
dashed line in (a) at different temperatures, (c) Current dis- 
tribution J y obtained by inversion of B tr ans image in (a), (d) 
Profiles of J y along the yellow line in (c) at different T with 
indicated irradiated and pristine regions. 

We now investigate the nature of the resulting liquid 
phase and search for a possible new transition between 
two different liquid phases. Since our DMO measure- 
ments do not show any FOT equilibrium magnetization 
steps within the liquid region of the phase diagram, we 
have developed a new differential imaging technique to 
search for a possible liquid-liquid transition. Instead of 
modulating the applied field or temperature, a transport 
current with periodically alternating polarity is applied 
to the sample and a corresponding differential image is 
acquired and integrated, typically over a hundred cycles. 
As a result, the self-induced magnetic field B trans gen- 
erated by the transport current according to the Biot- 
Savart law can be measured to a very high precision. 
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With this differential imaging we are able to determine 
the distribution of currents as low as 0.1 mA, about three 
orders of magnitude improvement compared to previous 
methods [18]. 

Figure 3a shows B trans in the central part of the sam- 
ple resulting from 50 mA transport current. This image 
displays the behavior at T = 74.2 K and B = 145 G, 
which is above the B^ D (T) line. Here the entire sam- 
ple is in the liquid phase and no field variation between 
the irradiated and pristine regions is observed by DMO 
using either field or temperature modulations. Surpris- 
ingly, despite the fact that both the irradiated and pris- 
tine regions are liquid, strong variations in B trans are 
found in Fig. 3b, revealing a highly inhomogeneous dis- 
tribution of transport current that reflects large differ- 
ences in the vortex dynamics. Such an ability to resolve 
spatial variations of vortex dynamics within the liquid 
phase is a unique feature of this new differential tech- 
nique. Using the matrix inversion [19,20] of B trans , a 
direct measure of the spatial distribution of the trans- 
port current J y is obtained as shown in Figs. 3c and 
3d. The current is strongly focused into the irradiated 
regions reaching densities up to seven times higher than 
in the unirradiated regions. This means that at this tem- 
perature the vortex liquid phase in the irradiated regions 
is different from the pristine homogeneous liquid; it has a 
much lower resistivity and the vortices experience a much 
higher effective viscosity. Figures 3b and 3d also show 
the evolutions of B tra ns and J y with increasing temper- 
ature. At Tdi = 79 K the current distribution becomes 
uniform quite abruptly and remains uniform at higher 
temperatures, showing that above Tdi the two liquids be- 
come identical. As discussed below, Tdi is thus the der- 
ealization temperature separating a nanoliquid phase, in 
which interconnected liquid droplets are caged by a solid 
skeleton, from a homogeneous liquid. Figure 2 shows the 
location of this derealization line Bdi(T) on the phase 
diagram, which is our central result here. 

We have also performed transport measurements si- 
multaneously with the DMO imaging. Figure 4 shows 
the measured four probe resistance R of the sample as a 
function of temperature at various fields. The arrows in- 
dicate the Tdi(B) as derived from the DMO self- field mea- 
surements. The resistance in the liquid phase of BSCCO 
has Arrhenius behavior [21,22] as seen in the inset of 
Fig. 4. We observe that in the homogenous liquid re- 
gion, R decreases smoothly with decreasing temperature 
until Tdi is reached. Below Tdi a nanoliquid is formed 
in the irradiated regions, significantly reducing R. In 
our special geometry, however, the volume of the irra- 
diated regions (brown regions in Fig. la) is less than 
half of the total volume of the sample. Therefore we 
should expect that, even if the resistance of the irradi- 
ated regions should become zero, the total measured R 
of the sample should drop approximately only by a factor 
of two, since the pristine regions remain liquid (blue in 



Fig. la) and maintain highly resistive Arrhenius behav- 
ior. Such a drop between two parallel Arrhenius lines is 
indeed observed and is emphasized by the dotted lines in 
the inset of Fig. 4. However, the striking observation is 
that the drop is by factor of more than 20. This behav- 
ior shows that the localization of the skeleton of vortices 
on the CDs significantly enhances the c-axis correlation 
as follows: Because of the high anisotropy of BSCCO 
the applied current flows in a very shallow layer under- 
neath the top surface of the crystal where the electrical 
contacts are attached [23-25] (see Fig. 4 lower inset). 
The characteristic depth of the current is inversely pro- 
portional to the effective anisotropy j p = (p c / pab) 1 ^ 2 , 
where p c and p a & are the c-axis and in-plane resistivi- 
ties. In BSCCO 7 P w 500 and may change considerably 
depending on oxygen doping and temperature. p c is de- 
termined by the degree of relative fluctuations and mis- 
alignment between the vortex pancakes in the adjacent 
Cu02 planes. The formation of the rigid skeleton aligns 
the pancake stacks residing on the CDs and also reduces 
the fluctuations of the vortices residing in the voids. As a 
result the c-axis correlation and the Josephson coupling 
between the pancakes in the nanoliquid phase should be 
significantly enhanced, giving rise to a large reduction 
in p c and anisotropy as compared to the homogeneous 
liquid. Consequently, the nanoliquid in the irradiated re- 
gions allows the current to penetrate much deeper into 
the crystal, thereby significantly increasing the effective 
thickness of the sample and reducing R. Since R drops 
by a factor of more than twenty, instead of a factor of 
two, we conclude that the nanoliquid is at least an order 
of magnitude more isotropic than the homogenous liquid. 
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FIG. 4. (color online). Sample resistance vs. T at various 
fields. The arrows indicate the derealization temperatures 
Tdi as derived from DMO self- field data. Top inset: example 
of Arrhenius plot of R(T) at B = 125 G with indicated Tdi and 
pristine melting . Bottom inset: image of the gold contacts 
on the crystal with irradiated regions shown schematically by 
the circles. 
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We now discuss the phase boundaries in Fig. 2 from 
the symmetry breaking standpoint. There are two sym- 
metries that are broken upon freezing into the vortex 
solid phase [1]: the formation of the Abrikosov lattice 
within the nanocrystals locally breaks the continuous 
translational symmetry of the liquid, and the appear- 
ance of the supercarrier density and the resulting abil- 
ity to sustain dissipationless currents along the vortices 
breaks the longitudinal gauge symmetry. At B^ D (T) the 
nanocrystals melt into nanodroplets, restoring the local 
transverse translational symmetry. The rigid skeleton, 
however, remains intact due to the strong pinning by the 
CDs, thus preserving the longitudinal c-axis correlation 
along the CDs up to the Bdi(T) line. The resulting nano- 
liquid phase, consisting of interconnected nanodroplets 
caged in a porous skeleton, has a unique property of ex- 
tended c-axis coherence together with the absence of in- 
plane critical current. 

Several recent theoretical studies and numerical simu- 
lations [7,10,11] indeed find evidence for superconducting 
coherence along the c-axis in such an intermediate liq- 
uid phase. More specifically, using the analogy between 
vortices and two-dimensional bosons [1,2], an intermedi- 
ate superfluid state in which the superfluid condensate 
of bosons (vortex liquid) coexists with localized bosons 
(vortices pinned on CDs) has been predicted [10]. With 
increasing T the intermediate vortex liquid system (nano- 
liquid) undergoes a sharp derealization at Bdi(T) above 
which all the vortices become delocalized (homogeneous 
liquid), thereby rendering the system incapable of sus- 
taining any supercurrents along the vortices. This re- 
sults in Bdi(T) oc exp(— T/To), where To describes the 
effective depinning energy scale for vortex derealization 
[10], estimated to be To ~ 5 K for typical parameters 
of BSCCO. The red lines in Fig. 2 are fits to the data 
resulting in To ~ 7.7 K for the B^ = 60 G sample and 
T « 13.9 K for the B^ = 20 G crystal (inset) which 
had a slightly higher T c . We have also studied a crys- 
tal with B^ = 47 G (not shown) giving T w 7.3 K. In 
all three cases the theoretical temperature dependence 
provides a very good description of the experimental 
Bdi(T) line. Note, however, that the theory does not take 
into account the discreteness associated with the pancake 
vortices in BSCCO. An alternative estimate, based on 
the Berezinskii-Kosterlitz-Thouless like scenario of der- 
ealization of two-dimensional pancake vortices [4], also 
yields an exponential dependence of Bdi(T) with proper 
parameters. Thus, we surmise that the localization of 
vortices on the porous skeleton depresses the c-axis re- 
sistivity and breaks the longitudinal gauge symmetry in 
the nanoliquid, which is restored only above the Bdi(T) 
line. It is commonly believed that the formation of the 
solid vortex lattice causes simultaneous breaking of the 
two symmetries across the FOT. Our findings here show 
the first phase diagram identifying two separate phase 
boundaries B^ D (T) and Bdi(T) across which the two 



symmetries of the vortex matter are apparently broken 
separately. This unique property is the result of the het- 
erogeneous nature of the vortex nanoliquid. 

It is interesting to note that at fields below the kink 
(marked by arrows in Fig. 2) the two symmetries are ap- 
parently broken simultaneously, where the B^(T) line 
closely follows the derealization line Bdi(T). At these 
fields the nanocrystals contain just a few vortices which 
remain solid as long as the caging skeleton exists. The de- 
localization of the skeleton at Bdi (T) thus causes instan- 
taneous melting of the entire lattice. This is a novel kind 
of vortex lattice melting in which derealization [4,10] of 
one sub-system of vortices induces a FOT in the entire 
heterogeneous system. Note that the finite distribution 
of the nanocrystal sizes results in broadening of the FOT, 
as indicated by the error bars in Fig. 2. The FOT appar- 
ently commences upon melting of the regions with larger 
crystallites while the global homogeneous liquid is estab- 
lished only upon melting of the smallest of the crystal- 
lites. Interestingly, this broadening is hardly resolvable 
for lower irradiation doses, as seen in the inset of Fig. 2. 

In conclusion, a low density of CDs is shown to signif- 
icantly perturb the vortex liquid resulting in formation 
of an intermediate nanoliquid phase in which liquid vor- 
tex droplets are caged by vortices localized on CDs. We 
find evidence for a new derealization line at which the 
nanoliquid transforms into a homogeneous vortex liquid 
forming a "Y" shaped phase diagram. The nanoliquid 
phase has unique properties resembling the long sought 
disentangled vortex liquid [26] with high c-axis correla- 
tion but zero transverse critical current. We conjecture 
that the nanoliquid phase with coexisting localized and 
delocalized particles may be generic to strongly interact- 
ing systems in the presence of disorder. 
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